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Abstract 

A linear theory of dust acoustic (DA) and dust ion-acoustic (DIA) waves in a self-gravitating 
collisional dusty plasma contaminated by variable-charge impurities is presented for a self-consistent 
closed system. The ion-drag force arising from the ion orbital motion as well as momentum transfer 
from all the ions and self-gravitation are taken into consideration. The physical processes, viz., dust- 
charge relaxation, ionization, recombination and collisional dissipations are also taken into account 
self-consistently. The generalized dispersion relation describing the coupling of the acoustic modes 
with the dust-charge relaxation mode (CRM) and damping of the waves is derived and analyzed 
numerically. 
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1. Introduction 



Charged impurities or dust grains frequently appear in space and laboratory plasmas [1,2]. 
Many low-temterature technological plasmas are contaminated by highly charged dust or impu- 
rities, since they carry a considerable ammount of negative charge of the plasma. The massive 
(compared to the ion mass) dust particles significantly affect the charge balance in plasmas, and 
producing the specific dust-charge fluctuations [3] lead to a new plasma mode, called charge relax- 
ation mode (CRM). There are many important factors which can significantly affect the dispersion 
characteristics of dusty plasma waves, for which we need to emphasize the open character of dusty 
plasma system demanding regorous treatment of collision mechanisms specific to dust-contaminated 
plasmas, esspecially those leading to plasma recombination. 

In the earlier investigations, it was shown that acoustic modes can be unstable in presence 
of ion-drag and ionization, and coupling between the DA and DIA modes is strong enough in 
the long- wave length run (see, e.g., Ref. 4). In the latter, DA and DIA modes interact through 
the dust-charge branch [5]. This is very different from the direct coupling between the long- wave 
length acoustic modes without the charge variation. However, in these investigations [4,5] the 
collisional effects are taken not in a self-consistent manner. Moreover, the ion-drag force which 
is of practical interest in processing plasmas are also treated as independent on dust charge. In 
absence of ionization, recombinatiion and collisional dissipations, the Langmuir waves were found 
to be unstable due to coupling with the CRM [6]. It was also found that the DAWs can be strongly 
affected by the particle creation and loss mechanisms [7]. The effect of dust-charge variation on 
Langmuir and ion-acoustic waves in a low-temperature plasma was investigated [8,9] by taking 
into consideration the self-consistent particle balance. It was shown that dust-charge relaxation 
is significantly affected by ionization and recombination. These processes maintain the averaged 
background particle number densities self-consistcntly during the perturbations by acting as sources 
and sinks. They also define the equilibrium or steady state. It was also shown that these dissipative 
processes due to inter-particle collisions as well as elastic Coulomb and the inelastic dust-charging 
collisions can lead to a net damping of the Langmuir waves [8] as well as waves in typical laboratory 
plasma [8,9]. 

It is known that the ion-drag forces are very important in the dispersion relation for modes in 
dusty plasma. These are discussed by a number of authors [10,11]. The ion-drag forces may arise 
from the ion orbital motion around negatively charged dust particles as well as from the momentum 
transfer from all the ions which are collected by the dust grains. The instability may arise due to 
the ion fluid compression by the ion-drag force. The latter, which is proportional to the square 
of the particle radius, dominates over the electrostatic force when the dust grains have sufficiently 
large sizes. 

The gravitational effects become important when the sizes of the dust grains become consid- 
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erable. For too small grains, trapping by magnetic field lines dominates much as it does for other 
plasma particles. Too large grains on the other hand follow gravitationally bound orbits, without 
being distracted from them by electromagnetic forces. Only for medium-sized grains both the ef- 
fects balance each other. In plasmas like protostellar clouds, there may be competition between the 
gravitational self-attraction and electrostatic repulsion between the grains. When self-gravitational 
interaction due to the heavier dust component is included, dusty plasmas are subject to macroscopic 
instabilities of the Jeans type. When a relative streaming between the plasma and dust component 
is present, Jeans-Buneman type instabilities are excited [12-14]. In the absence of equilibrium dust 
streaming and the dc electric field, the dust grains are held under the combined influence of the 
gravity force and the equilibrium drag force involving unperturbed ion flow towards the dust grain 
surface. 

In this report wc reconsider the problem of linear DA and DIA wave propagation in a self- 
gravitating dusty plasma with dust-charge variations, ionization, recombinations as well as colli- 
sional effects due to collisions between the electrons and dust, ion and dust etc. in a self consistent 
manner. We also consider the influence of the ion orbital drag force which depends on the ion 
velocity as well as dust charge. It is found that by these effects the DA and DIA wave propagations 
as well as instabilities are somewhat modified. 

2. The Model 

We consider the propagation of linear DA and DIA waves in a self-gravitating eoUisional dusty 
plasma with finite effective electron (Tg) and ion (Tj) teperatures. The size of the dust grains(rci) is 
assumed to be much smaller than the average intergrain spacing (d = (3/47rndo)^/"'^)), the electron 
Debye length (Xdc) and the wave length {rd << d << X£)e,k^^), so that they can be considered 
as massive point particles similar to multiply charged negative (or positive) ions in multispecies 
plasmas [15]. The charge of a dust grain may vary due to microscopic electrons and ion currents 
flowing into its surface, and these currents appear because of the potential difference between its 
surface and the adjacent plasma. We assume that the ion-drag force (Fcn) arising from the ion 
orbital motion around the dust grains dominates over that arising from the momentum transfer of 
all the ions which are collected by the dust grains. This is valid as long as the collection impact 
parameter is much smaller than the orbital impact parameter and the Coulomb logarithm far 
exceeds the unity. Thus the orbital motion related ion-drag force acting on the dust grains can be 
taken as [16] 

Fdi = 4:TmimiVitVibl/2HVit) (1) 
where Vn = (ff -|- SV^Jir)^/"^ is the total ion speed, 6^/2 = '^did^ I'^iVii is the orbital impact 
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parameter and A{Vit) = In (A^^g + ft^/2)/(^c + ^w) the Coulomb logarithm integrated over 

the interval from the collection impact parameter be = rd{l + 2Zdoe'^ /rdruiV^) to the electron Debye 

length. Other parameters will be discussed in the following. 

The dynamics of electrons, ions and negatively charged dust grains in presence of gravity and 

ion-drag force as well as ionization, recombination, attachments of electrons and ions to the dust, 

elastic and inelastic collisions and plasma particle collisions is governed by a set of equations: 

^ + V.{neVe) = -Vedrie + fion?le - (^effnl (2) 

I ,.e.f.f^^ I dne _ e d(t> 

dt e e m^n^ Qz rUe dz 

^ + V.{niVi) = -VidiH + v^ourie - Pefjnl (4) 

dt * ' miTii dz rrii dz 

^ + ^.{ndVd) = (6) 
^ Td dud Qd d(j) Fdi dip 

— + [Vd.VjVd + VdnVd H ^ = (7) 

dt nidTid dz rud dz rud dz 

= A-KGmdnd (8) 

V^(?!) = -AneiZiUi - rie + — n^) (9) 

e 

^ + -d% = -\le0\^ + \h0\^ (10) 

at Ueo riio 

where (l){ip) is the electrostatic (gravitational) potential and mj,Vj,nj (including the equilibrium 
value rijQ ) are the mass, velocity and density of the j-th species { j = i,e,d for ion, electron 
and dust ) respectively. Furthermore, G is the gravitational constant, Vd is the dust grain radius 
and Qd = —ZdS and ZjC are the charges of dust components and ions. In Eqs. (2) and (4) 
'^{e,i)d &re the attachment frequencies of electrons and ions to the dust grain so that —i^{e,i)d'^e,i 
represent their losses in densities; Vion is the rate of electron impact ionization of the neutral 
atoms, f3eff = (3 — Psi,f3 is the volume recombination rate and (3si is the stepwise ionization rate. 
In Eqs. (3) and (5) u^-^/ = Vei + ^{e,i)n + ^fi + effective collision frequencies. Here u^i is 

the rate of electron-ion collision, u^n is the a— particle (a = e, i) collision with neutral, uf,- are the 
rate of elastic (Coulomb) collisions with the dusts and arc the inelastic charging collision rates. 
Also in Eq. (7) I'dn is the effective frequency of the dust-neutral collision [17]. In the charging Eq. 
(10), uf^ = rdiOpiA/ \/2TTVTi is the charging rate of the dust particle [3], defined by the equilibrium 
electron and ion currents given by 
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lio = T^rjZiC Uioil ^ (12) 



where Acpg = qdo/C,C = rail + rd/Xoe) is effective grain capacitance and qdo the stationary- 



grain charge. Here A = I + t + Z,t = Ti/T^.Z = z^e^ /rdT^,VT(^e,i) = y {Te,i/me,i),ujp^i,d) = 
^ATrn(^i^d)o^id^'^ /''^i,d)- '^^^ rates of elastic electron and ion-dust Coulomb collisions are 

= rt^^^X^expiZ) = .t^j{A - l)exp{-Z) (13) 

where A is the Coulomb logarithm. The effective charging rates are [18] 

^ P{t + Z){A + Z) ^ 3 ;,n,o {t + Z){^ + Z) 

AZ 2 ' neo (2r + Z) ^ ' 

where P = ZdUda/rieo- Also the rates of electron and ion capture at the grain surface is [3,18] 

'^ed = ^^^^^l^t={^ + P>id (15) 

The expressions for the ionization, recombination rates as well as for electron-ion, ion-neutral col- 
lisions etc. can be had from Ref. [19]. 



3. Dispersion relation 

To obtain the self-consistent equilibrium state we assume that the pressure is not so low such 
that the recombination losses prevail over the diffusion losses, i.e., ambipolar diffusion can be 
neglected. In equilibrium, the system is quasi-neutral so that 

ZioUio = HeO + Zdorido (16) 
Prom Eqs. (2) and (4) we obtain the stationary electron and ion densities as 

«eO = Z ' '^iO = neO (17) 

Peff m 

Wc note that the ionization degree must be high enough so that fion > i^ed^ otherwise no self- 
consistent stationary state exists. 

Assuming the perturbations to vary as exp(ikz — icut) and linearizing the Eqs. from (2)-(10) 
we obtain for the perturbed electron density and velocity {nj = njo + n^i, j = e, i; Zd = Zdo + Zdi) 

iek'^cj) ^ ^l^ed \ ^ 



riei 



(18) 



= ^. eff ' (19) 
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where rje = lo — i{yed — ^2^ion) — ^^Tel i'^ + i^l^^) ) ^2 = 2 — ^1 . Here depends on the used model 
of the direct ionization, = corresponds to density-fluctuation independent ionization, whereas 
^1 = 1 that for density-fluctuation dependent. 

For the ion fluid velocity and perturbed ion density we have 



Vi = 



where 



Ai = 



+ 



eff' 



mi{u + rieme{LO + ii'e^-^) 



eff 
i 

c = 

A2 = 



(20) 



(21) 



riio rie 



dVid \ _ f dved \ 
dZdJza, \e\dzj 



\Ieo\k^(pA3 



For the dust charge and number density perturbation we have 

Zdi = -i 

ridok 



ridi 



Tnid 



rie{uj + ivf-) 



(22) 
(23) 



Here 



A.. 



+ 



irjeZi 



1 + iC/Vi 
rUe {ijj + iiyl^^) r]imi (to + iiy^-^-^ ) 



ill 



eOl 



er]e 



\dZdJz,, m \dZdJz, 



(24) 



Also r = + j^jy^^ _ k'^V^^ - ujjd,V^^ = Td/md,bid ^ vrZ^oeVSTj, Aq = A{Vito) and uj'^j = 
AirGmdndo is the Jeans' frequency. Finally, substitutions of Eqs. (18)-(23) into Eq. (9) lead to the 
dispersion relation for DA and DIA waves 



OJ. 



pe 



I • eff 



pi 



eff 



Ci 



(25) 



where 



Ce = - 

Ve 



1 + ^ — (6i^j< 

Vi 



^l^ed) 



Ci = k'^/Vi 



_ 2 , V^^dOJpd ■ 47re|/eo|fc=^A3 



fc2 [dZd)z,o ' Vi \dZd). 



ndo 



v = — 

rrid 



o.^ Ni/2 VTibf.Ao eZi k^V^i / . 
qao + i8{2Tr)^/^niomi ''^J — + ^ eA 



Vi 



|Je0|A2A3 



(26) 
(27) 



Equation (25) describes the linear coupling of Langmuir, DA and DIA waves with the CRM u) = 
—ii^d^ where the eff'ects of the variation of electron and ion capture rate by the dusts are taken 
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into consideration. Since it is very difficult to analyze Eq. (25) eitfier analytically or numerically 
keeping all the terms in order, we rewrite the equation in a comparatively simplified form on the 
assumption that 

« ^{e,i)d (28) 



However, for larger dust grain size and higher ion number density, the condition (28) can be violated. 
Note that one can recover the dispersion relation equation of Ref. [8] for Langmuir plasma waves. 
Here we only discuss the DA and DIA modes. Thus, neglecting the the variations of v^^. ^^^ with 
respect to we obtain from (25) 



(29) 



Here 



D{u;,k) = l---f 



^2. 



,2 



A,, 



A, 



X={^'' + irio; + Ci - fcV|jA„ K = {u' + ir2u; + C2 - k'V^i) 



2t/2 



A. = ^ 



<^ld^id 



CO + Zfj 



eff 



1 + 



X 



{lo^ + iT^uj + C4 - kW-f^ 



Cl = K^H'^^ed - ^2l^ion), C2 = mi'f^ 



(3 = (2 + Z.Cll^f^M", C4 = C2 + — 

^ed l/j ■' ■' 



r4 



2z^ed) 



ri + — {^2l^ion - 'ifed) , D^h (w) = W + iuf 



me Zil/gd 
V 

-I 

X 

Ts = r4 H -ZiFi, 

rrii 



Pcoupii^, k) = + ^rgo; + Cs - k^V^i - ^Ziev^e)col, 



rrii 

S5 = (.4 ACl) ^ = 

rrii erieo 



Note that u is the correction of the dust charge relaxation frequency arising from the perturbation 
of the electron density. 

Equation (29) describes the coupling of the acoustic waves given by D(a;, A;) =0 with the CRM 
given by D^hioj) = 0. We discuss the two modes separately as below by considering the hierarchy of 
the frequency scales for the DIA and DA respectively as Icjd/aI ~ 10^ — and \ujda\ ~ 1 — 10^/s. 

I. Dust ion-acoustic modes 



Assuming kVn, kVrd « lo « kVre and Td « Te,Ti we obtain from Eq. (29) the following 
dispersion relation for DIA modes 



iff 



+ CU^J + iuj{T2k'V4, + TsCU^e) + (C2 - ^pi)k'V^e + CsCU, 



Dch+ 
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ik^V^^uld^id{J^ + iV2U + C2)l Dch + iu?ul^uDiff{u? + iV^u + Cs - —Zik^V^^) = (30) 

where Dijf = {u + iu^^^). Setting lo = loi + id'i + id", where ooi = oOpikXoe, we find from (30) in 
the hmit k'^Xjjg << 1 

r/ QR1 + PR2 , . 

^' = - P^ + Q-^ (^'^ 

„ QR2 - PRi 
= - P^ + Q^ ^''^ 

where 

P = bujfai — Sufa^ + 05, Q = 6a;f — Auifa^ + 2a;ia4 

D 6 4 I 2 D 5 3 I 

Ri = uji — ujia2 + ujia4 — ae, R2 = aiUi — ujia^ + ujia^ 
ai = -(^5 + + r^f^^ + Ta), as = -Cs + + z^^^^) + (^^^^ + ^s) + ^^^^^rg + k'Xl^^i 
as = HC5 - i^t^r,) + ufiCs - Tsi^t") + Csi^t" - k'Xh {^>f^ + ^t) - ^la^rd) 
a4 = -{K^ + ^fCs>f^ - k'Xl, [i^tiu^au^id - oj^uf^) + uld'^iaT2 

05 = k^X])^ujlaUJid{ufT2 - C2), 06 = X]y^ul^ujidvf C2 

In Eq. (31) negative sign indicates that a frequency downshift takes place for the case considered, 
and comparing the values of all dissipative terms one can show from Eq. (32) that in presence 
of ion-drag force there may be an increase in the damping decrement. In a dusty argon plasma 
with Tg 10 eV, Ti 1 eV, T„ ~ 0.015 eV, nio/n,,^ ~ 10, 5fim, Z^o ~ 10^ we find for 

^2 = 2, 6[ « -2.393 X 10^, 6'{ « -1.54 x 10^ in absence of ion-drag and 6[ ^ -2.3545 x lO^s, S'{ « 
-1.5737 X lOVs 

in presence of ion-drag. That is the ion-drag force leads to decreasing of the 
amount of shift and the growth rate of damping. From the detailed numerical analysis we also find 
that the amount of shift increases with the electron-to-ion temperature as well as ion-to-electron 
density ratio. 

The Eq. (28) is higher order in the wave eigen frequency, for experimental practice one can at 
once find the spatial attenuations by setting kXoe = Ki + iK2 as 



^2 ±{XdXm + YdYn) + ^{Xl + Yl) {Xl + Y^) 



Here 



Xb = a;'xi - iyfx2 + Auf^, Xm = u^u^^ - vf^Xi) 
Yn = u{vfxi + X2-^), YN = u\uffxs + XA), A = ^Z.^a;^^^^ 

Tfli 

2/ 2 I A 2 \ 2-n e// 2 -n 

= u;^r2 + co^uf^ico^ + C2 - a;^J + u^ld^^idioj'' + C2) 

= ^2(^2 ^ _ ^ch^^ _ -^5), ^4 = j,ch^^2 ^ ^ ~(^2 ^ 
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For a typical dusty plasma the DIA wave frequency is ~ 10^ /s, so that the self-gravitational effect 
on DIA modes being smaller can be neglected. We have plotted Ki and K2 vs real u for different 
values of electron to ion temperature and ion to electron density ratio (Figs.(l), (2)). We find 
that both the ratios increase the damping decrement of the wave. But when the density ratio 
takes the value 6.0 the propagation stops in a short range of frequency domain and the maximum 
decay rate is then 4.6 x 10~^ (Fig. 2 dashed line ). Further increase of ion density increases the 
damping decrement. We numerically verify the condition |A;^A|)g| << 1 and the restriction u < Upi 
for DIA modes. To estimate the ion number density we find that for a typical argon plasma pa- 
rameters described above the threshold ion number density (n^Q*^^*) is ~ 9.0 x 10^^ cm~^ when 
u 1.07 X 10^/s,Ki Ri 0.291, if2 ~ 0.953. A slight decrease of n^o = 8.95 x lO^^cm"^ would lead 
to a increase of Ki and K2 {Ki ^ 0.293, K2 ~ 0.958) and hence a violation of the condition for kXoe- 

II. Dust acoustic modes 

Assuming kVrd « uj « kVre, KVn we obtain from Eq. (29) the following dispersion relation 
for DA modes 



D. 



ch 



+u{u + iuf^){ir5U + Cs - k^V^i - "^Zik^V^jDjd = (34) 
where 

Djd = [UJ + lUJVdn - ^jd)^ V=^+ T ^ 

^i^ed 

A detailed numerical analysis shows that frequency downshift occurs due to coupling of DA modes 
with the CRM. 

From Eq. (34) we obtain for kXce = + iKj as 



Kr 



(35) 



Ki « -^^JYdYn (36) 



Here 



w^T^Ai - A2A3] + ulaUJid {oj'^ii^f^ - riujid) + yfrni^f^Ti) 
Yd = uldUJid {^id^f^V - ^fi^f^ - Wid)) + c^pi(??A2 + A1A3) 
Xn = w^AiA4 - A2A5, Yn = A1A5 + A4A2 

A 2 2 cff 

V = ^ + r^ , Ai =a; -w-d- 

A2 = Uj'^iydn + l^f^i'^'^ - C^ld)' A3 = ufr) + V{1 + Uid/ZiUed) 
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A4 = Ca - yfr^ - ^^rs, As = Czi^f + + u?V^ 

Note that in the dispersion relation (34) we have included the self-gravitation due to the heavier 
dust components, and in the derivation of Eqs. (35) and (36) we assumed that X]:,,n » Yd,n- 
For an argon plasma described earlier we find that u^-^/ ~ 10^/s,fd„ ~ 10~^/s, fjon so 
that the above approximation can be made. 

Prom Eq. (34) we find that Xd > or < according as 

This shows that both the gravity and ion-drag force stabilize the instability of the DA modes. 
The numerical calculations for (35) and (36) are carried out for the parameters of a dusty argon 
plasma and plotted against real oo as shown in the Figs. (3)- (5). To show clearly the influence 
of gravity on the DA modes we numerically calculate an estimate for n^o in presence and ab- 
sence of self-gravity to verify the restrictions uj < Upd and |A;^A^g| << 1. We find that for the 
parameters nio/ueo = 10.0, Tg/Tj = 10.0, r,^ ~ 5/Ltm,n„o ~ 3.0 x lO^^cm"^ the threshold for n^o 
is ^ 13.92 X lO^cm-3 when oj ^ 0.335/s,ifi ^ 1.83 x 10'^, K2 -5.89 x 10"^ below which the 
DA modes do not exist in presence of gravity. Whereas in absence of gravity the upper limit for 
Udo for the existence of DA modes is 13.05 x lO^cm"^. We also find that in presence of gravity 
as the dust number density increases from the threshold value, the instability of the modes tend 
to decrease (Fig. 4, the solid, dotted, dashed and dadottcd lines), e.g., at a particular frequency 
cj = 0.24/s, the values of K2 are -5.48 x 10-^,-4.88 x 10"^-3.79 x 10"^ and -1.22 x 10"^ for 
Hdo = 13.9 X 10'' cm"^, 15.7 x lO''' cm"^, 17.4 x lO'^ cm~^ and 19.1 x 10'' cm~^ respectively. That is, 
gravity influences to stabilize the instability in the dense region of dust grains. Figure (3) shows that 
for a particular dust density where lOj^ = 7.29 x lO^^s^^, if we increase the value of Tg/Ti from 1.2 
the value of Ki{> 0) also increases, i.e., electron and ion temperature also play an important role to 
the damping decrements of the modes in presence of gravity. We have also shown that for particular 
values of temperatures and densities (Tg/Tj = 1.6, Udo = 8.7 x lO"^ cm~^ , Uio / UgQ = 10.0), Kj again 
becomes negative but in absence of gravity and its magnitude increases with the frequency (Fig. 5) . 



Conclusions 



Wc have presented a general model for linear acoustic waves as well as Langmuir waves sclf- 
consistently with the presence of dust grain charge variation as well as the influence of ion-drag and 
self-gravitation. The model also includes the ionization, volume recombination, electron and ion 
capture by the dusts, as well as Coulomb elastic and inelastic charging collisions of the particles. 
From our general dispersion relation one can recover Langmuir wave as in Ref. [8], DA and DIA 
waves as in Ref. [4,5] in some particular considerations and also the DA waves as in Ref. [7] for 
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stationary neutrals. We have shown that in presence of ionization, recombination and all other 
collisions, linear coupling of both DA and DIA waves with the dust charge mode (CRM) leads to 
damping and a frequency downshift of the waves. This is typical consequence of the coexistence of 
acoustic modes and CRM in dust-contaminated plasmas [9] . The presence of ion-drag force which 
is opposite to the electrostatic force decreases the amount of shift and the damping increment of 
the waves. We have presented graphically the dependence of real and imaginary part of DA and 
DIA wave number on the wave eigen frequency for a dusty argon plasma for different values of 
temperatures and densities. We have shown that the instability of the acoustic modes are modified 
somewhat by the influence of ion-drag and self-gravity force. Prom our numerical analysis we have 
estimated the ion and dust number densities due to the restrictions u) < Upi and uj < Upd for DIA 
and DA modes, and found that in the region of dense dust grains the gravity influences to stabilize 
the instability of the DA modes. Due to our lack of experimental data we could not compare our 
results with those of existing experiments. Our results will be more applicable to many laboratory 
as well as astrophysical plasmas. 
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Figure Captions 



Figure (1): Plot of the real (a) and the imaginary (b) parts of kXoe vs the DIA wave frequency 
u ~ (Eq.(30)) in an argon plasma. The solid (0 < a; < 118.72), dottted (0 < a; < 119.44), 

dashed {0 < u < 116.96) and dadotted (0 < a; < 118.6) lines correspond to Te/Ti = 1.0,3.13,5.0 
and 19.2 respectively. The other parameter values are given in the text. The corresponding fre- 
quency intervals indicating those where the waves exist. For uj = 0,Ki =0 and K2 ~ 0. 

Figure (2): Plot of the real (a) and the imaginary (b) parts of kXoe vs the DIA wave frequency 
u ~ lO^/s (Eq.(30)) in an argon plasma. The solid (0 < w < 77.0), dottted (0 < a; < 164.8), 
dashed (0 < a; < 2.8) and dadotted (0 < a; < 236.0) lines correspond to riio/neo = 1.0,4.0,6.0 and 
8.0 respectively. The other parameter values are given in the text. The corresponding frequency 
intervals indicating those where the waves exist. As u ^ 0, Ki, K2 0. 

Figure (3): Plot of the real (a) and the imaginary (b) parts of /cA/jg vs the DA wave frequency 
cj ~ 1/s (Eq.(34)) in an argon plasma. The solid (0 < a; < 0.21), dottted (0 < a; < 0.23) and dashed 
(0 < a; < 0.24) lines correspond to T^/Ti = 1.2, 1.6 and 2.0 respectively Here 00]^ = 7.29 x lO^V^^ 
and the other parameter values are given in the text. The corresponding frequency intervals indi- 
cating those where the waves exist. For a; = 0, i^i « 10~^ and K2 = 0. 

Figure (4): The influence of gravity on DA modes: 

Plot of the real (a) and the imaginary (b) parts of kXue vs the DA wave frequency a; ~ 1/s (Eq.(34)) 
in an argon plasma. The solid {0 < lo < 0.335), dottted (0 < u; < 0.332), dashed {0 < lo < 0.314) 
and dadotted {0 < lo < 0.249) lines correspond to lo^^ = 0.1166,0.1312,0.1458 and 0.1604/s2 re- 
spectively. The two lower solid lines of both (a) and (b) are corresponding to no self-gravitation. 
The other parameter values are given in the text. The corresponding frequency intervals indicating 
those where the waves exist. For co = 0,Ki and K2 ~ 0. 

Figure (5): The same as the dotted line in Fig. (3). Here Ki becomes negative in absence of 
gravity, showing that instability increases with the frequency. 
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This figure "Fig. 1. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/nlin/0501016vl 



This figure "Fig. 2. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/nlin/0501016vl 



This figure "Fig. 3. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/nlin/0501016vl 



This figure "Fig.4.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/nlin/0501016vl 



This figure "Fig. 5. GIF" is available in "GIF" format from: 



http://arXiv.org/ps/nlin/0501016vl 



